Density functional theory study of the electronic structure of fluorite Cu 2 Se 
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We have investigated the electronic structure of fluorite Cu2Se using density functional theory 
calculations within the LDA, PBE and AM05 approximations as well as with the non-local hybrid 
PBEO and HSE approximations. We find that Cu2Se is a zero gap semiconductor when using either 
a local or semi-local density functional approximation while the PBEO functional opens up a gap. 
For the HSE approximation, we find that the presence of a gap depends on the range separation 
for the non-local exchange. For the occupied part in the density of states we find that the LDA, 
PBE, AM05, PBEO and HSE agree as regards the overall electronic structure. However, the hybrid 
functionals result in peaks shifted towards lower energy compared to the LDA, PBE and AM05. 
The valence bands obtained using the hybrid functionals are in good agreement with experimental 
valence band spectra. We also find that the PBE, PBEO and HSE approximations give similar 
results regarding bulk properties, such as lattice constants and bulk modulus. In addition, we have 
investigated the localization of the Cu d-states and its effect on the band gap in the material using 
the LDA+C7 approach. We find that a sufficiently high U indeed opens up a gap but that this U is 
unphysically large. 



I. INTRODUCTION 

The copper chalcogenides CU2X (X = S, Se, and Te) 
are of possible technological interest because of their 
thermo- and photoelectric properties as well as ionic con- 
ductivity. Cu 2 _ a: Se, which is the focus of this study, has 
in particular received attention due to a high ioni c co n- 
ductivity, with possible applications in solar cells^l as 
well as a good material for thermoelectric convertersP 
At room temperature Cu2_zSe is a rather good p-type 
conductor with an optical band gap of 1.23 eVP The 
Cu2-a;Se system has a rather complicated atomic struc- 
ture where the phase diagram consist of two phases: The 
low temperature a-phase and the high temperature (3- 
phaseP The /3-phase with the space group Fm3m has the 
Se atoms in a face-centered cubic (fee) environment while 
the superionic Cu atoms are randomly distributed on in- 
terstitial positions in the structure, where the majority of 
the Cu ion s ar e positioned close to the tetrahedral inter- 
stitial sites!^ The low temperature phase is stable up to 
about 400 K and has a lower symmetry crystal structure, 
where the Cu atoms are not diffusing. We point out that 
the a to (3 transition varies with the stoichiometry and 
for x = 0.15 — 0.25 the high temperature phase is stable 
at room temperatureP In addition, we note that both 
the high and the low temperature phases can be imag- 
ined as deviations from the fluorite crystal structure with 
the difference lying in the occupation of the tetrahedral 
interstitial positions. 

Recently Liu et alP showed that Cu2- x Se exhibited 
excellent thermoelectric properties for a bulk material 
with a thermoelectric figure of merit zT = 1.5 at 1000 K. 
The reasons for the favorable thermoelectric properties 
appear to be a low thermal conductivity due to the 
quasi-liquid behavior of the superionic Cu atoms com- 



bined with a rather high value of the power factor S 2 a, 
where S is the Seebeck coefficient and a is the electrical 
conductivity!^ It is interesting that a material with a very 
simple chemical formula has such a low thermal conduc- 
tivity since this is usually reported for complex systems 
containing heavy elements, for example in doped skut- 
tcrudites and clathrates, or in nanostructures. 9 Further 
investigations of the physical properties of Cu2-a;Se as 
well as other similar systems are therefore important in 
order to obtain a fundamental understanding of the ther- 
moelectric properties of the Cu2_a;Se system in particu- 
lar and superionic solids in general. In this paper we will 
present results obtained by density functional theory cal- 
culations on the electronic structure of Cu2Se. The aim 
with the present study is to obtain a theoretical under- 
standing of the electronic structure of Cu2Se by means 
of density functional theory (DFT) calculations in order 
to establish the accuracy and applicability of various ap- 
proximations for the exchange-correlation energy (XC) 
functional within DFT. 

The paper is arranged as follows: In Section [TT] we will 
present the necessary details of our calculations, in Sec- 
tion III we will present our results and in Section IV we 
will summarize our results and draw conclusions. 



II. DETAILS OF THE CALCULATIONS 

We have performed density functional theory calcu- 
lations for Cu2Se in the fluorite structure, space group 
Fm3m, shown in Fig. [TJ This structure is an ideal- 
ized version of the high temperature phase of Cri2Se 
and serves as a good model for determining the accuracy 
of various approximations for the XC functional. As is 
shown in Fig. [T] the fluorite structure has the Se atoms 
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FIG. 1. The crystal structure of fluorite Cu2Se where Cu 
atoms are depicted by small (white) spheres and Se by large 
(black) spheres. The Cu-Se bonds are shown in grey. 



positioned at a regular fee crystal lattice while the Cu 
atoms occupy the tetrahedral interstitial positions in the 
fee lattice. 

The Kohn-Sham equation has been solved using the 
projector augmented wave methoc^as it is implem ented 
in the Vienna ab initio simulation package (VASP)f 13 * 14 ! 
Calculations have been performed using the local den- 
sity approximation (LDA)J^the generalized gradient ap- 
proximation of Perdew, Burke and Ernzerhof (PBE) 16 
as well as with the AM05 functional of Armiento and 
MattssonP^H^ Since local and semi-local approximations 
to the XC functional are known to perform badly when 
determining the band gaps of semiconductors and insula- 
tors we have in addition performed calculations with the 
more advanced hybrid functionals of Perdew, Ernzerhof 
and Burke (PBE0)pSl an d of Heyd, Scuseria and Ernzer- 
hof (HSE) J2H2H These functionals have been found to be 
an improvement of the LDA and PBE when regarding 
for example bond distances and dissociation energies in 
molecules as well as structures and band gaps of bulk 
materials EiO Interestingly, the AM05 has been found 
to yield the same level of accuracy as the hybrid func- 
tionals when regarding bulk properties of crystals, such 
as lattice constants and bulk modulus.^ 

The PBEO and HSE are non-local hybrid theories 
where the exchange part of the exchange-correlation 
functional includes some non-local contribution. In the 
PBEO approximation the XC energy can be expressed by 
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where E PBE and E RBE are the exchange and correlation 
terms in the PBE approximation and E BF is the Hartree- 
Fock exchange energy. In the HSE approximation the XC 



where LR and SR denote long range and short range 
parts respectively, a is a mixing parameter governing the 
strength of the non-local exchange and uj is a screening 
parameter that controls the spatial range over which the 
non-local exchange part is important. The correlation 
energy in the HSE approximation, E BSE , is taken to be 
identical to the PBE correlation energy as in the PBEO 
approximation. The amount of Hartree-Fock exchange 
in the HSE is the same as in the PBEO, i.e. a = 1/4, 
while the range separation parameter, uj, has to be de- 
termined by comparison with experimental data. It has 
been found that uj = 0.2 — 0.3 A" 1 give good results re- 
garding structural as well as electronic properties of ma- 
terials, with uj = 0.2 A -1 being the optimum choice .^H^ 
Unless otherwise specified, uj = 0.2 A -1 has been used in 
the HSE calculations. We note that for uj = Eq. ^ is 
equivalent to the PBEO and, furthermore, Eq. g asymp- 
totically reaches the PBE for to — > oopD We will return 
to this point in Sec. |III| 

An issue with the use of common density functional 
approximations is the description of strongly localized 
orbitals such as d- and f-orbitals. For this reason we have 
also investigated how the localization of the Cu 3d-states 
affects the electronic structure as well as the band gap in 
the system by employing the use of an on-site Coulomb 
interaction within the LDA+[7 approximation, using the 
approach of Dudarev et alP^We note that in the Dudarev 
et al. approach the on-site Coulomb interaction, U, and 
on-site exchange interaction, J, do not enter separately 
but only in the combination of an " effective U" , U e t / = 
U — J. However, throughout the presentation we will 
simply refer to the "effective U" as U. The U has been 
applied solely to the Cu 3d-states. 

For the main part of the calculations, the plane wave 
basis set was cut-off at 1000 eV and we have used a k- 
point mesh of 30x30x3(P3for the primitive cell of Cu 2 Se. 
However, since the hybrid functionals are much more 
computationally cumbersome than the local and semi- 
local approximations, we have used a plane wave cut-off 
of 800 eV and a k-mesh of 10 x 10 x 10 in combination with 
the hybrid functionals. 

The pseudopotential (PP) for Cu has the 3d and 4s 
electron states treated in the valence while the semi-core 
3s and 3p states have been placed in the core of the PP. 
For Se the PP has the 4s and 4p electron states treated 
as valence states. 
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FIG. 2. (Color online) Calculated density of states and 
projected density of states obtained by using the PBE. The 
upper panel is the total DOS and the DOS projected on the 
Cu and Se atoms. The middle panels is the projected Cu 
DOS with Cu 3d- and 4s-states marked. The lower panel is 
the projected Se DOS with Se 4p- and 4s-states marked. The 
vertical dashed line marks the position of the Fermi level Ep. 
In the middle panel the DOS for the Cu 4s-states has been 
magnified 10 times. 



III. RESULTS 

The lattice constants and bulk modulus obtained us- 
ing the various functionals are summarized in Table [I] 
We note that the PBE, PBEO and HSE functionals give 
very similar results in good agreement with experimental 
observations, whereas the LDA overbinds. The AM05 is 
the functional that gives the lattice constant closest to 
the experiment with a deviation that is less than 1 %. 
When adding the on-site Coulomb interaction, U, to the 
LDA we find that the lattice constants decreases with in- 
creased U. The bulk modulus also becomes more similar 
to the PBE, PBEO and HSE results as U increases. Usu- 
ally, the application of an on-site interaction is expected 
to increase the lattice constant since the U reduces the 
kinetic energy of the electrons by allowing them to lower 
their energy by localization. This leads to a lessened 
binding energy for the electrons for which a U has been 
applied and in turn an increase in the lattice constant. 
A reduced lattice constant may nevertheless appear since 
the population of strongly bonding s- or p-states may in- 
crease at the expense of the d-state occupation, resulting 
in an increased binding strength. 

In Fig.[2j we show the electron density of states (DOS) 
as evaluated using the PBE functional. It is clear from 



TABLE I. Calculated bulk properties of fiuorite Cu2Se ob- 
tained using the LDA, PBE, AM05, PBEO and HSE approxi- 
mations for the XC functional, a is the lattice constant and B 
is the bulk modulus. The lattice constants and bulk modulus 
are evaluated using a Birch-Murnaghan 3rd order equation of 
stateP^The lattice constants and bulk modulus in the second 
and third columns correspond to the XC functional in the first 
column while the fifth and sixth columns correspond to the 
U given in the fourth column within the LDA+t/ approxima- 
tion. 

XC functional a (A) B (GPa) U (eV) a (A) B (GPa) 
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the results in Fig. [2j that the valence band consists of 
three regions, a lower region between 14 eV and 15 eV 
below the Fermi level, which solely is built up by Se 4s- 
states, a middle region in-between 5 eV and 7.5 eV below 
the Fermi level that is a mixture of states derived from 
both Cu and Se, and an upper region from the Fermi 
level down to 3.5 eV that mostly consists of Cu 3d-states 
even though there are mixing with the Se 4p-states. The 
Cu 4s-states are also present in all three regions. We 
note that there is a Cu 3d peak in the DOS at about 
3 eV below the Fermi level that is positioned in-between 
two regions of Cu 3d and Se 4p hybridized regions, the 
lower region between about -8 eV to -5 eV and the up- 
per region from about -2 eV to the Fermi level. This 
behavior of a Cu 3d peak centered in-between regions of 
hybridized Cu 3d and Se 4p states has been found exper- 
imentally by Domashevskaya et al.pSl even though the 
exact positions of the peaks differs which we will turn to 
shortly. We also note that the result presented in Fig. [2] 
is in agreement with the previous tight binding calcu- 
lation of Garba and Jacobs^ as well as with the LDA 
result of Kashida et alPJ A notable result is that the 
PBE functional yield a zero gap between the occupied 
valence bands and the unoccupied conduction bands. As 
was mentioned in Section [TJ optical studies have reported 
a band gap of 1.23 eVP This is a rather small gap and 
density functional theory is known to underestimate the 
band gap of semiconductors and insulators when using 
local and semi-local approximations. In order to investi- 
gate the band gap for this system further we will have to 
resort to more complex approximations to the exchange- 
correlation functional. 

In Fig.[3]we compare the DOS obtained using the LDA, 
PBE, AM05, PBEO and HSE as well as LDA+C/. We 
find that the overall shape of the DOS is very similar 
to the PBE, which was discussed previously, except for 
when very large [/-values have been applied. In the fol- 
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FIG. 3. (Color online) Calculated density of states for PBE, AM05, PBEO and HSE in the left column and for LDA and 
LDA+U, with U — 4 eV, 8 eV and 10 eV, in the right column. The vertical dashed lines marks the position of the Fermi level 
Ef for each case. The total DOS is shown in black and we also show the projected DOS on the two different atomic species; 
Cu in red and Se in blue. 



lowing discussion we divide the valence states into lower 
(Se 4s-states), middle (mainly Cu 3d- and Se 4p-states) 
and upper (including the main Cu 3d peak) valence re- 
gions. Compared to the PBE the predominant trend is a 
shift of the states towards lower energies. For the LDA 
and AM05 this shift is rather small. For the PBEO and 
HSE we find that the lower valence region is found be- 
tween 15 eV and 16 eV below the Fermi level, compared 
to 14 eV to 15 eV for the PBE. The middle valence re- 
gion is also shifted downwards with the same amount 
(~ 1 eV). For the upper valence region the lower edge is 
shifted downwards with about 1 eV for the hybrid func- 
tionals. For very large U the upper valence states are 
shifted more than the states in the middle valence region 
and for U — 8 eV and 10 eV the upper and middle va- 
lence regions from the PBE description have merged into 
a single region. We note that the valence band DOS ob- 
tained by the PBEO and the HSE functionals are in very 
good agre ement with measured x-ray photoelectron spec- 
tra (XPS) ! 11 ! 29 ! The main peak in the spectra at 3.5 eV 
come from Cu 3d-states and is excellently reproduced by 
the main Cu peak obtained in the DOS from the PBEO 
and the HSE approximations as well as for LDA+C7 with 
U = 4 eV, see Fig. [3} For larger [/-values the Cu d-states 



are found at too low energies. The LDA, PBE and AM05 
have the Cu d-states at slightly to high energy compared 
to the experiment. As regards the conduction bands and 
especially the band gap, we find that all approximations 
have a region of rather small DOS above the Fermi level. 

Figure [4] shows the calculated band structure for the 
upper valence band as well as for a few conduction bands 
along high symmetry directions in the Brillouin zone. As 
is clear, only the PBEO and LDA+U with U = 10 eV cal- 
culations reveal a band gap. In all other cases the highest 
occupied molecular orbital (HOMO) and the lowest un- 
occupied molecular orbital (LUMO) touch at the T-point. 
The LDA, PBE and AM05, has an almost flat conduction 
band along the A line (r to L), where the LUMO at L is 
only ~ 0.1 eV above the Fermi level. For these function- 
als the lowest conduction band are degenerate with two 
valence bands at T. The same is also true for the HSE 
approximation as well as for LDA+U with [7 = 4 eV and 
8 eV. However, as can be seen in Fig. |4j the HOMO-2 
band at the T-point behaves differently depending on the 
XC functional. This band is for the PBE band structure 
found at 1.19 eV below the Fermi level, while for the HSE 
calculation it is found at 0.07 eV below the Fermi level. 
For the PBEO this band has moved up to the Fermi level 
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FIG. 4. (Color online) Calculated band structures along directions of high symmetry in the Brillouin zone for Cu2Se obtained 
using the LDA, PBE, AM05, PBEO and HSE functionals. The vertical solid lines label high symmetry coordinates while the 
horizontal dashed line marks the position of the Fermi level, Ep. Note the common feature that in order for a gap to be open 
the HOMO-2 band, marked by the (red) dashed band, has to reach the Fermi level, which only happen in the case of the PBEO 
and LDA+I7 with (7=10 eV. The fat bands show the amount of Cu 4s character of the lowest conduction band and the three 
highest valence bands. 



and the HOMO at the T-point is triply degenerate with 
a direct gap of 0.47 eV. The same situation is the case 
for U = 10 eV but the gap is 0.30 eV. 

When an on-site Coulomb exchange is added to the 
Cu 3d-states, the main effect is that the Cu 3d-states are 
pushed further down below the Fermi level as witnessed 
by the projected DOS in Fig. [3] We find that the upper 
valence bands in Fig. [4] is becoming broader and the con- 
duction bands are shifted towards higher energies. The 
major effect in the band structure of Fig. [4] is a shift in 
the low- lying conduction band along the A-line (r to L) . 
In order for a gap to be open very large [/-values are 
required. The lower right panel of Fig. [4] shows the re- 
sulting band structure with U — 10 eV. The behavior 
of the valence and conduction bands with U = 10 eV 
is rather similar to the PBEO results in the Fermi level 
near region, i.e. about -2 eV to 2 eV, with a direct band 
gap at r of similar size as the PBEO. We find that the 
HOMO-2 band is key for the formation of a band gap 
within the LDA+C/ approximation. At U = 8 eV, this 
band just reaches the Fermi level and marks a transition 
situation and a small increase of the Coulomb exchange 



interaction to 8.1 eV (not shown) gives a small gap at T. 

Common features of the band structures shown in 
Fig. [4] are that the highest valence band, as well as the 
lowest conduction band, are found along the A line and 
that the ordering of the bands is similar irrespective of 
approximation used. We also note that the valence bands 
obtained by the PBEO and HSE are very similar, as is 
also witnessed by the DOS in Fig. [3] The difference be- 
tween these two hybrid approximations lies in the ob- 
tained conduction bands, with PBEO these bands are 
found higher above the Fermi level compared to the HSE. 
The band structure data shown in Fig. [4] are quantified 
in Table [H] where we show the evaluated energy differ- 
ences between HOMO and LUMO levels at various high 
symmetry points in the Brillouin zone. 

We have shown that the existence of a band gap is 
closely related to the position of the HOMO-2 level at 
the r-point: If this band reaches the Fermi level a gap is 
opened. Common for all calculations shown in Fig. [4] is 
that the triply degenerate state at the Fermi level at T is 
derived from hybridization between Cu 3d-states with t 2g 
symmetry and Se 4p-states. The HOMO-2 band at T is 
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TABLE II. Calculated transitions between HOMO and LUMO levels at high symmetry points in the Brillouin zone. Energies 
are given in eV. 
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derived from hybridized Cu 4s- and Se 4s-states, however, 
if a gap is open this is no longer the case: Whenever a gap 
is open the HOMO-2 band reaches the Fermi level and 
the Cu 4s and Se 4s character is found in the LUMO. This 
can be seen in Fig. [4] where we show the amount of Cu 
4s character of each of the bands close to the Fermi level. 
As shown in Fig. [4j there is a region close to T where the 
HOMO-2 band has a rather large Cu 4s character in the 
PBE, AM05 and LDA band structures. In the HSE band 
structure the Cu 4s character of this band is much more 
narrow about the T-point, however, it is still present at T. 
For the PBE0 bands there is no Cu 4s character for this 
band at T. The same transition of the band character 
is found for the LDA+U calculations, where U = 8 eV 
is a transition state with Cu 4s character remaining in 
the HOMO-2 at V which is completely transferred to the 
LUMO for U =10 eV. 
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FIG. 5. (Color online) Evaluated band gaps using various 
values for the range separation parameter oj within the HSE 
approximation. Note that tj — in the figure is the PBE0 
result. Furthermore, as ui — > oo the PBE result is obtained. 
The points for u = 0.01, 0.05, 0.1, and 0.3 A -1 have been eval- 
uated using the lattice constant obtained with ui = 0.2 A -1 
and the points are connected by a fitted spline. The vertical 
dashed line marks ui = 0. 



Since the PBE0 and HSE calculations yield very sim- 
ilar valence bands but differ significantly regarding the 
positions of the conduction bands, it is relevant to an- 
alyze the differences between these two approximations. 
In Fig. [5] we show the calculated band gaps obtained by 
varying the range separation parameter, cj, within the 
HSE approximation. The standard values for this pa- 
rameter are 0.2 — 0.3 A -1 , where w = 0.2 A -1 is con- 
sidered to be the optimum choicel^iHUl We find that by 
going below the optimized u window, a gap opens up, and 
the PBE0 value is recovered at w = 0. The reason for 
this behavior can be linked to that in the HSE approxi- 
mation, long range exchange contributions are neglected, 
see Ref. [5TJ By writing the PBE0 expression in Eq. 
in short and l ong range parts along the lines of the HSE 
approximation-- 1 we find that 
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where a is the amount of Hartree-Fock exchange included 
(1/4 in the present case). Thus the difference between 
PBE0 and HSE is the range of the included long-range 
exchange. By varying u> it is possible to bring the HSE 
result in closer agreement with the PBE0, see Fig. [5| 



IV. DISCUSSION AND CONCLUSIONS 

We have performed density functional theory calcula- 
tions of the electronic structure of fluorite Cu2Se using a 
variety of approximations for the XC energy functional. 
We find that the PBE, PBE0 and HSE approximations 
give very similar results regarding lattice constants and 
bulk modulus. The hybrid functionals PBE0 and HSE 
both show excellent agreement with experimental valence 
band spectra, although all functionals give rather simi- 
lar results regarding the overall electronic structure, see 
Fig. [3] In addition, we find that the PBE0 is the only 
functional that yield a non zero band gap. It is also 
possible to open a gap by adding an on-site Coulomb 
interaction to the Cu 3d-states. However, we conclude 
that localization of the Cu 3d-states cannot account for 
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the formation of a gap since this requires unphysically 
large [/-values. Furthermore, compared to experimental 
valence band spectra such large [/-values push the Cu 
3d-states too low. 

We also find that the presence of a band gap is dictated 
by the relative position of the Cu 4s and Se 4s hybridized 
band and the Cu 3d and Se 4p hybridized band at the 
T-point. When a gap is open, i.e. for the PBEO and for 
very large U, the Cu 4s and Se 4s hybridized band is 
found in the LUMO above the triply degenerate HOMO 
that consists of Cu 3d and Se 4p hybridized bands. In 
all other cases the Cu 4s and Se 4s hybridized band is 
positioned below the Fermi level. 

In comparison to experimental band gaps we conclude 
that the present calculations do not reproduce the ob- 
tained optical gap of 1.23 eV obtained by Sorokin et alP 
However, this can be due to several reasons where one 
is related to the structure of the system. It is possible 
that defects, most likely Cu vacancies or interstials, af- 
fects the electronic structure in such a way that larger 
gaps are obtained and further investigations along these 
lines are under way. We point out that the experiment 
was performed at room temperature and for the ordered 
low temperature phase of Cu2Se which has a lower sym- 
metry than the fluorite structure. For Ag2Te, which is 
similar to Cu 2 Se, GGA calculations give a zero gap for 
the fluorite structure, while a gap opens up for the mon- 
oclinic structure.^ We speculate that the same may also 
be true for Cu2Se, i.e. that a sizeable gap opens for the 
low temperature phase and that the fluorite structure has 
a vanishingly small gap. 

Furthermore, our calculations reveal that the band 



structure around the T-point is such that it would not 
create a clear peak in the optical spectrum but rather an 
energy independent low intensity signal which can easily 
be confused with a background. In addition, our DOS 
calculations shown in Fig. [3] reveal a region with very 
small DOS above the Fermi level irrespective of the XC 
energy functional that has been used. 

In addition, it has been found that HSE is better com- 
pared to the PBEO at describing the band gap in small 
gap semiconductors, where the latter is found to over- 
estimate the size of the band gapP^ Thus, in summary, 
despite the fact that PBEO is the only functional opening 
up a gap, we believe that HSE with uj ~ 0.2 A -1 is the 
functional that more correctly describes the electronic 
structure of Cu2Se. In order to further analyze the opti- 
cal spectra by theory it will be necessary to resort to solv- 
ing the Bethe-Salpeter equation or using time-dependent 
DFT where optical gaps can be described properly due 
to the inclusion of the electron-hole interaction in these 
methods. 
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